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SUMMARY 


Interior noise in the fuselage of a twin-engine, propeller-driven aircraft 
with two propellers rotating in opposite directions is studied analytically. 

The fuselage is modeled as a stiffened cylindrical shell with simply supported 
ends, and the effect of stringers and frames is averaged over the shell surface. 
An approximate mathematical model of the propeller noise excitation is formu- 
lated which includes some of the propeller noise characteristics such as sweep- 
ing pressure waves around the sidewalls due to propeller rotation and the local- 
ized nature of the excitation with the highest levels near the propeller plane. 
Results are presented in the foirm of noise reduction, which is the difference 
between the levels of external and interior noise. The influence of propeller 
noise characteristics on the noise reduction is studied. The results indicate 
that the sweep velocity of the excitation around the fuselage sidewalls is 
critical to noise reduction. 


INTRODUCTION 

One of the main sources of interior noise in propeller-driven aircraft is 
the propeller noise excitation. A main path by which noise is transmitted to 
the interior is through the fuselage structure. The present paper is devoted to 
the study of the noise reduction provided by the fuselage structure when sub- 
jected to propeller noise excitation and the influence of propeller noise char- 
acteristics on the noise reduction. 

Some information is available on the nature of the propeller noise and its 
near field behavior (refs. 1 to 7) which reveals a few characteristics of pro- 
peller noise such as (1) the pressure field rotating with the blades causing a 
sweeping action around the sidewalls of the fuselage, (2) the localized nature 
of the noise field on the fuselage surface with highest levels near the pro- 
peller plane, and (3) the noise spectrum dominated by propeller harmonics. Ref- 
erences 8, 9, and 10 analyze interior noise in fuselages using statistical 
energy analysis and deterministic analyses for exterior noise pressure inputs 
such as turbulent boundary layers, jet noise, or plain acoustic waves. Interior 
noise in propeller-driven light aircraft with flat fuselage sidewalls has been 
studied experimentally (ref. 11) , and analytically considering uniform pressure 
to represent propeller noise input on the sidewalls (ref. 12) . 

The present investigation extends previous analyses by employing a more 
realistic representation of propeller excitation pressures and complements the 
studies discussed in reference 12 with regard to fuselage geometry by consider- 
ing cylindrical fuselages. An approximate mathematical model has been developed, 
using information available in the references, for the propeller noise excita- 
tion which includes some of the propeller noise characteristics. Twin-engine 
aircraft with propellers rotating in opposite directions are considered. The 
noise fields from the propellers on either side of the fuselage are assumed to 
be symmetrical about the central vertical plane of the fuselage. The objective 



of the present investigation is to study the quantitative influence of the pro- 
peller noise characteristics on the noise reduction through the fuselage side- 
walls* Two sizes of fuselage structure are studied, one representative of a 
large passenger aircraft fuselage and the other corresponding to that of a light 
aircraft fuselage. 

The overall development of the analysis in this study closely follows that 
in reference 10. The fuselage is modeled as a cylindrical shell stiffened by 
stringers and ring frames and having simply supported ends. The effects of 
stiffeners are averaged over the fuselage surface. The analysis for the interior 
noise in the fuselage is carried out using the approximate model of propeller 
noise input on the fuselage structure and the normal mode approach to determine 
the response of the fuselage and the resulting pressure fluctuations in the 
interior acoustic field. The frequency range considered is from 0 to 500 Hz and 
includes the blade passage frequency and the first few harmonics, which are the 
dominant sources of propeller noise. The noise reduction through the fuselage 
sidewalls is the difference between the excitation noise level on the fuselage 
surface and the space-averaged interior noise level within the fuselage. Results 
showing the influence of propeller noise characteristics on noise reduction are 
presented. 


SYMBOLS 

The units used for physical quantities defined in this paper are given in 
the International System (SI) of units. 

Ar^As cross-sectional areas of ring frame and stringer, respectively, m^ 

aj^,a 2 fb 2 /b 2 distances defining the fuselage surface area subjected to pro- 

peller noise excitation, shown in figure 1, m 

constants appearing in equation (B2) 

C . modal constant defined in equation (Bll) 

mni 

c speed of sound in air, m/sec 

E modulus of elasticity of fuselage material. Pa 

E^,Eg moduli of elasticity of ring frame and stringer materials, 

respectively. Pa 

^mni generalized force for (mni) mode 

f frequency, Hz 

fj. reference frequency, Hz 

Gr#Gs shear moduli of ring frame and stringer materials, respectively. Pa 

H frequency response defined in equation (6) 
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h thickness of fuselage material, m 

^mni integrals defined in equations (11) , (12) , and (13) 

modified Bessel function of first kind and of nth order 
moments of inertia for ring frame and stringer, respectively 
Bessel function of first kind and nth order 
Jj.,Jg torsional constants for ring frame and stringer, respectively 

k wave number 0 )/c defined in equation (B3) 

L length of fuselage, m 

lirs differential operators in equation (Al) , r,s = l, 2, 3 

^r'^s ring frame and stringer spacings, respectively, m 

generalized mass of (mni) mode 
moment resultant in x-direction 

m, 5 number of half waves along cylinder axis 

stress resultant in x-direction 

n, r) number of full waves around cylinder circumference 

P pressure amplitude of excitation pressure. Pa 

p interior noise pressure. Pa 

Pg excitation pressure. Pa 

dynamic magnification factor for acoustic modes 

Qmni dynamic magnification factor for (mni) structural mode 

^mni generalized coordinate of (mni) structural mode 

R radius of fuselage, m 

r radial coordinate of fuselage, m 

/\ 

r dimensionless radial coordinate r/R 

s number of radial acoustic waves 

t time, sec 
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a 


mni 


u 


Ui,U2 


V 

mnx 


Vi,V2 


W 


mni 


w 




z # z 
r' s 


a 

"^mni 

0 

li 

P 

Pa 

Po 


Pr'^s 


$ . 
mni 




mode functions for fuselage displacements in x-direction 

fuselage displacement in x-direction, m 

fuselage displacement amplitudes in x-direction, m 

sweep velocity of exciting pressure waves, m/sec 

mode functions for fuselage displacements in y-direction 

fuselage displacement in y-direction, m 

fuselage displacement amplitudes in y-direction, m 

mode functions for fuselage displacements in z-direction 

radial displacement of fuselage, m 

radial displacement of air particle in fuselage interior, m 
coordinate in longitudinal direction of fuselage, m 
coordinate in circumferential direction of fuselage, m 
coordinate in radial direction of fuselage, m 

eccentricities of ring frame and stringer centroids from fuselage skin 
middle plane, m 

parameter defined in equation (B3) 
quantity defined in equation (B17) 

phase angle of frequency response of fuselage displacement, rad 
Poisson's ratio for material of fuselage 
density of fuselage skin material, kg/m^ 

average density of skin-stringer fuselage structure, kg/m^ 
density of air, kg/m 

densities of ring frame and stringer materials, respectively, kg/m^ 
mode functions of (mni) acoustic modes 

coordinate in circumferential direction for interior acoustic 
space y/R 

quantity defined in equations (B4) 
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0 ) 


frequency, rad/ sec 


03^ structural modal frequencies, rad/sec 

cOinns acoustic modal frequencies, rad/sec 

Subscripts : 

i, j ,m,n,r, s,r)f 5 integers 

Dots over symbols denote time derivatives. 

Primes over symbols denote derivatives with respect to spatial coordinates. 
Bars over symbols denote amplitudes. 


ANALYSIS 

Since propeller noise is predominant near the blade passage frequency and 
the first few harmonics, which lie in the low-frequency range, a fuselage model 
that is simple but suitable in the low-frequency range is required. Such* a 
model is available in reference 13 in which the effect of stringers and ring 
frames are averaged over the fuselage surface. Shown in figure 1(a) is a 
fuselage structure idealized as a simple cylindrical shell with reinforcing 
stringers and ring frames. The fuselage is excited by propeller noise excita- 
tion, and the resulting vibrations of the fuselage structure excite the interior 
air space causing interior noise. The displacement response of the fuselage 
walls is calculated using a normal mode approach. The interior noise in the 
fuselage is obtained by equating the displacement of the fuselage wall to that 
of the air particle near the wall. This plan of development follows refer- 
ence 10 closely. 

The following discussion gives the main features of the analysis. Details 
of the derivation of the forced-response equations for the fuselage structure 
are given in appendix A, and details of the interior noise derivations are 
given in appendix B. 


Displacement Response of Fuselage Structure 

The equation of motion for the (mni) mode of the fuselage structure is 
(see eq. (A13) ) 


•• ^ . 2 /4-\ - ^mni^^^ 

^mni ^ Q ^mni ^ ^ ^mn^mni ^ ^ 


■^mni 




( 1 ) 


where represents the generalized coordinates and represents the 

natural frequencies of the fuselage structure. The subscripts (mni) specify 
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the resonant deflection shapes of the cylindrical shell where m=l, 2 , 3 , . . 
is the number of elastic half waves along the cylinder axis, n=l, 2, 3, . . , 
is the number of elastic full waves around the cylinder circumference, and 
i = 1 , 2 indicates the orthogonal solutions for mode shapes to specify the 
complete response since there is no preferred location of the origin of the 
y-axis. The first few modes are illustrated in figure 2. The dynamic magnifi- 


cation factors at resonance are denoted by Q„ 


The term 


is the gen- 


eralized mass of the (mni) structure! mode. (See eq. (All).) The term 
is the generalized force of the (mni) structural mode given as 


i(t) = Pe(x,y,t) 


where P 0 (x,y,t) is the external pressure loading on the fuselage surface and 
Wjnni(x,y) denotes the (mni) structural mode function. The complete radial 
displacement response of the shell is obtained by summing up the response of 
all the normal modes of the shell (eq. (a7c) ) as 






Assuming that the generalized force is harmonic at frequency o); that is 


^mni^^^ ^mni® 


the solution of equation ( 1 ) can be written as 


_ H (03/(0^) 

^mni^^^ 2 

^mn^\ini 


exp [-0(03/0)^)] 


where 


- 




9(u/u^) - tan-1 ) ^ ^ 





I 


Subsequently, the radial displacement of the shell is obtained by combining 
equation (5) and equation (3) as 

H((jO/u) ) j- — _ . . 

w(x,y,t) = > 


Propeller Noise Excitation 

Previous studies on the nature of the propeller noise (refs. 1 to 7) indi- 
cate that: 

1. The noise field generated by the propeller blades rotates with the 
propeller. 

2. The propeller noise exciting the fuselage surface has the highest 
levels near the propeller plane. 

3. The spectrum of the propeller noise is dominated by propeller harmonics. 

By employing these characteristics of propeller noise excitation an approximate 
mathematical model is formulated. 


In figure 1(b) the location of the propellers is shown relative to the 
fuselage. The propeller noise on the fuselage is illustrated by the hatched 
area around the fuselage. It is assumed that the propellers on either side of 
the fuselage are rotating in opposite directions and also that their noise 
fields incident on the fuselage sidewalls are symmetrical with respect to the 
fuselage vertical plane of symmetry. It is also assumed that the noise field 
is traveling around the fuselage with constant velocity V, as shown in fig- 
ure 1(b), by the arc around the fuselage. This approximate representation of 
the propeller noise field is mathematically expressed, for a particular pro- 
peller harmonic, as 


Pe (x,y,t) 


P exp 


ico 1 1 



Pe(.x,y,t) 


= P exp 


10 ) 





Pe(X/y/t) = 0 


(elsewhere) 
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where O) is the frequency of the harmonic, P is the amplitude of the exciting 
pressure, V is the velocity of the pressure waves traveling around the fuse- 
lage due to the propeller rotation, and a^^, a 2 /b 2 ^,b 2 specify the fuselage sur- 
face area subjected to propeller noise excitation. (See figs. 1(a) and 1(b).) 
The sweep velocity V and the frequency 03 can be related to the propeller 
parameters. In equations (9) , v = corresponds to the case of uniform pres- 
siare loading on the fuselage surface. 

Combining equations (9) and equation (2) and making use of equations (A8) 
and (A9) , the generalized force for the (mni) mode is obtained as 


F . (t) = F .e = Pe (I .) 
mni mni mni 


( 10 ) 


The integrals Imni equation (10) are evaluated to obtain 


I 


mni 


0 


( 11 ) 


due to the symmetry of the forcing pressure field about the vertical plane of 
symmetry of the fuselage, and 


^mn2 = (mira^/L) - cos (m7ra2/L)] 



( 12 ) 
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when n/R ^ 00 / V and 


^mn2 = (m^aj^/L) - cos (m7Ta2/L)] 


r 

sin 


nb^ 


^(b2 - bi)J 


nbi 


sxn 


R V 


R , ,, , > “^1 

+ (b2 - bj^) cos 


+ 1 




r 

cos 


nbi 

L— * v“=2 - ‘=1>J - ^ 


R V 


+ (b2 - b^) 


sin 



(13) 


when n/R = o)/V. In equation (13) when n/R = 0)/V, the sweep velocity V 
matches with the bending wave velocity around the fuselage 0)R/n and hence 
corresponds to the case of coincidence. Substitution of the generalized force 
given in equation (10) into equation (8) yields the displacement response of 
the fuselage structure . 


Interior Noise Analysis 

The internal acoustic field, which is excited by the fuselage vibrations, 
is defined by the wave equation 




(14) 


and the boundary conditions chosen for the ends and wall of the cylindrical 
section. To permit simplified analysis, open-end conditions are assumed; that 
is, the acoustic impedance, which is defined as the ratio of acoustic pressure 
to particle velocity, is taken as zero. While it is recognized that in an 
actual fuselage the ends have finite impedance, the use of the open-end condi- 
tions should have little effect on the calculated modal density of the sound 
field in the frequency range considered. 

An expression for the interior noise pressure is obtained by assimiing the 
boundary conditions at the wall to be such that the acoustic and the structural 
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analyses are coupled. The radial deflections of the shell are equated to the 
radial displacements of the air particles adjacent to the shell walls in the 
interior space. (See appendix B for detailed analysis.) The expressions for 
the space -averaged me an- square interior noise is given by (see eqs, (B16) , 
(B18), and (B25) ) 



In equation (15), (a) 0 is the condition of acoustic resonance. The 

acoustic resonance frequencies are given by (see eq. (B14) ) where 

m = 1, 2, 3, . . . is the number of longitudinal half waves, n = 1, 2, 3, . . . 
is the number of circumferential full waves, and s=l, 2, 3, ... is the 
number of radial waves. Nodal patterns for a few acoustic modes are shown in 
figure 3. At resonance the interior noise is given by 



(16) 


RESULTS 

Results are presented for two fuselage structures whose dimensions are 
shown in table I. Fuselage 1 is representative of a large passenger-aircraft 
fuselage and fuselage 2 corresponds to a light-aircraft fuselage. The prop- 
erties of the fuselage material and air in the fuselage interior are given in 
table II. The magnification factor for structural modes Qj^± is chosen as 50, 
which corresponds to 1 percent of the critical damping. This is the nominal 
value of damping for many actual structures. The magnification factor for 
acoustic modes is chosen as 25, which corresponds to 2 percent of the 

critical damping. Reference 14 suggests this value of the magnification factor 
for acoustically treated enclosed spaces in the low-frequency range. 

In order to provide a partial check of the structural part of the present 
analysis, natural frequencies of a simply supported, externally stiffened shell 
were calculated and compared with experimental results from reference 15. The 
experimental model had longitudinal stiffeners integrally built into the shell. 
As shown in figure 4, the agreement between analytical and experimental fre- 
quencies for this check case is good. 
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Natural Frequencies of Structure and Acoustic Field 


Analytical natxiral frequencies of both fuselages are shown in figure 5. 

In calculating the results, all the modes in the frequency range studied are. 

included (0 to 500 Hz) ; however, to simplify figure 5, not all of the modal 

frequencies are included. Figure 5 shows that the structural modal spectrum 
is significantly different for the two aircraft. For example, the lowest nat- 
ural frequency is about 28 Hz for fuselage 1 and about 108 Hz for fuselage 2, 

and fuselage 1 has more modes (approximately 250) than does fuselage 2 (13) in 

the frequency range. 

The acoustic mode frequencies of the fuselage interior space are shown in 
figure 6. As in the case of structural modes, all the acoustic modes in the 
frequency range studied were included in the calculations (0 to 500 Hz) ; how- 
ever, not all modes are shown in figure 6. Fuselage 1 had more acoustic modes 
in the frequency range studied (1200) than fuselage 2 (47) . The fuselage 
interior space receives acoustic energy in its acoustic modes, which are excited 
by the fuselage structural modes. For the particular end conditions assumed, 
the mnth acoustic mode will respond only to a structural mode having the same 
m and n numbers. (For other end conditions, the coupling between structural 
and acoustic modes is more complicated, as discussed in ref. 10.) In addition, 
when the natural frequencies of such mode pairs are close together, noise 
transmission into the fuselage interior takes place very efficiently. From the 
complicated nature of the modal frequency curves shown in figures 5 and 6, it 
can be seen that many such frequency matchings may occur. A search of the tabu- 
lated values of the modal frequencies revealed that approximate frequency 
matchings occurred for fuselage 1 at the following values: 


m 

n 

s 


^mns 

1 

1 

0 

46.3 

51.3 

1 

6 

0 

218.7 

204.1 

2 

6 

0 

218.0 

205.0 

2 

8 

1 

388.6 

384.0 

3 

2 

0 

88.0 

89.7 

3 

8 

1 

388.0 

385.0 

4 

8 

1 

387.5 

386.2 

5 

8 

1 

387.5 

387.7 


For fuselage 2, such frequency matchings occurred only at 


m 

n 

s 

mn 

^mns 

2 

2 

0 

269.7 

270.2 


The relation of these frequency matchings to the noise reduction will be 
discussed later. 
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Noise Reduction 


Numerical results for the interior noise are calculated from equations (15) 
and (16) and are presented in the form of noise reduction, which is defined as 


Noise reduction = 10 



(17) 


2 . 2 
where is the mean-square excitation pressure and p is the space-averaged 

mean-square interior noise pressure. The frequency range considered in the 
present study is from 0 to 500 Hz, since this range covers the blade passage 
frequency and the first few harmonics which are the predominant frequencies 
exciting the structure. In the results shown in figures 7 to 12, noise reduc- 
tion is studied as a continuous function of excitation frequency, for a number 
of sweep velocities. It is noted that an interior noise spectrum actually con- 
sists of tones at the blade passage frequency and its harmonics. The noise 
reduction at those corresponding frequencies can be obtained from the continu- 
ous noise reduction plot. Sweep velocity is a function of the propeller rota- 
tional speed and the distance from the propeller center to the fuselage side- 
wall. Hence, when the propeller rotation or the propeller center relative to 
the fuselage sidewalls are altered, the sweep velocity is affected and a differ- 
ent noise reduction plot must be computed. 

To illustrate the effect of blade passage frequency on noise reduction, a 
plot of noise reduction against frequency is shown in figure 7. The blade 
passage frequency is denoted as fj_ Hz. The noise reductions at and its 

higher harmonics are marked with circular symbols. If the blade passage fre- 
quency is changed to fj (either by changing the number of blades in the pro- 
peller or the speed of rotation of the propeller) but the sweep velocity V is 
kept constant, the resulting discrete frequency spectrum of noise reduction 
associated with f^^ and its harmonics changes as indicated by the square 
symbols in figure 7. Sweep velocity can be held constant by adjusting the 
distance between the propeller center and the fuselage sidewalls. It is seen 
from figure 7 that, by changing the blade passage frequency, the noise reduc- 
tion changes significantly, depending on whether it and its harmonics lie near 
resonances or away from resonances. Hence, noise reduction is very sensitive 
to blade passage frequency and structural and acoustic resonances. 

The effect of sweep velocity on the noise reduction is shown in fig- 
ures 8(a) and 8(b) for fuselage 1 and fuselage 2, respectively. The entire 
surface area of the fuselages is exposed to noise excitation as shown in fig- 
ure 8. Noise reduction is plotted against nondimensional frequency, with the 
reference frequency f^^^ as the nondimensionalizing frequency. The reference 
frequency is the frequency of m = 1, n = 0, shell mode obtained as described 
in appendix A. This reference frequency is approximately equal to the ring 
frequency discussed in reference 14. It is seen from figure 8(a) that, when 
the excitation is uniform (sweep velocity v = ^ in eqs. (9)), noise reduction 
is a minimum at the reference frequency, whereas when the exciting pressures 
are sweeping around the fuselage sidewalls the noise reduction at this frequency 
is considerably higher. At frequencies below f/fj- = 0.5, however, noise 
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reduction is higher for uniform excitation than for the case of excitation with 
sweep velocity. For the case of uniform excitation, covering the full circum- 
ference of the fuselage (that is, b 2 = Hr and = 0) , the integrals Ij^ni 

are zero except when n = 0, and hence, for uniform excitation, only n = 0 
structural modes are excited and transmit noise to the fuselage interior. The 
dip in the noise-reduction curve at the reference frequency for the case of 
uniform excitation in figure 8(a) may be explained by the large number of 
structural modes near this frequency as compared to other frequencies, as shown 
in figure 5. The dips in this curve below the reference frequency are due to 
the presence of acoustic modes only. In figure 8(b), for the frequency range 
shown, the noise reduction for uniform excitation is always higher than it 
is when sweep velocity is present. A comparison with figure 8(a) shows that 
this behavior is similar to that of fuselage 1 in the frequency range where 
(f/f^) < 0.5. From figures 8(a) and 8(b), it is seen that the effect of sweep 
velocity on noise reduction is quite significant, often making a difference of 
25 dB in comparison with uniform excitation. 

The variation of noise reduction with sweep velocity is shown in fig- 
ures 9(a) and 9(b) for fuselage 1 and fuselage 2, respectively. Only a portion 
of the fuselage surface area is exposed to excitation, in the vicinity of the 
propeller plane, as shown in figures 9(a) and 9(b). It is seen from these. fig- 
ures that, when the sweep velocity is reduced from 341.4 m/sec to 136.6 m/sec, 
the noise reduction increases by about 10 dB. These velocities are chosen 
due to the fact that many propellers have tip speeds in this range. The same 
trend of noise reduction, decreasing with increasing sweep velocities, has been 
observed to be generally maintained over a broad range of sweep velocities. 

This trend is shown in figures 10(a) and 10(b), for fuselage 1 and fuselage 2, 
respectively, where the noise reduction is plotted against sweep velocity at a 
few representative resonance frequencies. 

The effect of changing the axial length exposed to noise excitation on 
noise reduction is shown in figures 11(a) and 11(b) for fuselage 1 and fuse- 
lage 2, respectively. It is seen from both the figures that at most frequen- 
cies the noise reduction is greater by as much as 10 dB when a smaller portion 
of fuselage length is subjected to excitation. An exception is shown in fig- 
ure 11(b) at about 270 Hz where the noise reduction is substantially less when 
a smaller length of fuselage is excited. As indicated previously, the struc- 
tural and acoustic modes with m = 2 and n = 2 have nearly the same fre- 
quencies, resulting in a matching of their resonances and consequently a low 
noise reduction. This low noise reduction does not occur when the whole fuse- 
lage is excited because the m = 2 axial modes are orthogonal to the uniform 
excitation distribution. It is also interesting to note that the dips in the 
noise-reduction curve in figure 11(a) correspond with the modal frequencies 
discussed earlier where structural and acoustic modal frequencies match. 

The effect of circumferential pressure distribution on noise reduction is 
shown in figures 12(a) and 12(b) for fuselage 1 and fuselage 2, respectively. 

In these figures noise reductions are plotted for two different circumferential 
areas exposed to noise excitation, keeping the axial portion exposed to noise 
and the sweep velocity constant. A comparison of the two curves in figure 12(a) 
shows that the noise reduction does not vary substantially with circumferential 
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area for most frequencies* In figure 12(b) noise reduction is higher by as 
much as 20 dB when the circumferential area exposed to noise is reduced. These 
results indicate that, in general, it is important to have a knowledge of the 
distribution of noise in order to determine properly the noise reduction in a 
fuselage. 


CONCLUDING REMARKS 

The effect on noise reduction of propeller noise characteristics has been 
studied for a twin-engine, propeller-driven aircraft with two propellers rotat- 
ing in opposite directions. The fuselage is modeled as a cylindrical shell with 
stringers and ring frames, whose effects are averaged over the fuselage surface. 
The mathematical model of the propeller noise excitation used in this study 
includes some of the propeller noise characteristics such as pressure waves 
sweeping around the fuselage sidewalls due to propeller rotation and the local- 
ized nature of excitation with the highest levels near the propeller plane. 

The influence of the propeller noise characteristics on the noise reduction 
is studied. 

The results of this investigation show that sxibstantial changes in noise 
reduction occurred when sweep velocity and area of fuselage subjected to noise 
excitation were changed. Noise reduction generally decreased as the sweep 
velocity increased. Reducing the area exposed to excitation generally increased 
noise reduction. For the range of variables studied, differences in noise 
reduction of the order of 10 to 20 dB were observed. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 13, 1978 
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APPENDIX A 


EQUATION OF MOTION OF THE FUSELAGE STRUCTURE 

In the differential equation of dynamic equilibrium for the fuselage 
structure, Donnell-type strain displacement relations for the cylinder (ref. 13) 
and beam- type strain displacement relations for the stiffeners are considered. 
Inplane inertias are neglected since they do not significantly influence the 
radial modes of the shell, which are of interest in this investigation. The 
equations of motion for an element of the smeared shell can be written as 


Lj^lU + ^ 




13^ 

2 s'' 

33 W 


0 

0 


p;,(l - 


9^w 


Pe(5C,Y/t) 

Eh 


(Ala) 

(Alb) 


(Ale) 


where u, v, and w are shell displacements in the x-, y-, and z-directions , 
respectively, h is the shell thickness, E is the modulus of elasticity of 
the shell material, y is Poisson's ratio, t is the time, and Pg(x,y,t) is 
the excitation pressure on the fuselage surface. The term is the average 

density of the fuselage material given by 


Pa 


P + 


Ps^s 

h 


p^A 


r""r 


h I 


(A2) 


The differential operators in equations (Al) are as follows (see ref. 13) 


^11 


1 + 


E^Asd - y2) 
EhZg 


8^ 

9x^ 



(A3a) 


^12 ■ ^21 


/ I + V\ 3^ 

\ 2 / 8x 3y 


(A3b) 


^13 “ ^31 


}j 9 ^s^s^s(Z- ~ 9^ 

R^‘ Ehlg 9 jj3 


(A3c) 
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92 

^22 - 2 J 

-H 

'ax2 

^ EhZ^ 

^ 

9 y 2 

T — T ~ 

E A (1 - y 2 ) 

1 J- 

19 

ZrErArd - U^) 

^23 ~ ^32 

JL r 

EhZr ] 

R 3 y 

EhZ^ 


(A3d) 


(A3e) 


= ^ V4 + i 

TO’’ ' O 


^33 - 12 




1 + 


ErArd - U^) 


EhZ, 


Esd - U^)(ls + ZgAg) 34 ^ E^d - y^)(lr + ^r^r) 9^ 

EhZ„ ' 


9x^ 


2z^E^A^d - y2) J^2 

9y2 


EhZ, 


,2 s /G„J„ G„J„ 


3y" 


^ ^ (1 - I -s^s . 


REhZ, 


Eh 


3x^9y^ 


(A3f) 


With simple supports at each end of the shell, the boundary conditions to be 
satisfied at x = 0 and L are 


w = M^ = v = Nx = 0 


(A4) 


where is the moment resultant and is the stress resultant in the 

x-direction. The boundary conditions in equation (A4) are satisfied by 


u(x,y,t) = ui(t) 

COS 

(m7Tx/L) 

sin 

(ny/R) 

(A5a) 

v(x,y,t) = vi(t) 

sin 

(mTTx/L) 

cos 

(ny/R) 

(A5b) 

w(x,y,t) = wi(t) 

sin 

(mTTx/L) 

sin 

(ny/R) 

(A5c) 

by 






u(x,y,t) = U 2 (t) 

cos 

(mTTx/L) 

cos 

(ny/R) 

(A6a) 

v(x,y,t) = V 2 (t) 

sin 

(m7Tx/L) 

sin 

(ny/R) 

(A6b) 
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w(x,y,t) = W 2 (t) sin (m7Tx/L) cos (ny/R) (A6c) 

where m is the number of longitudinal half waves, and n is the number of 
circumferential waves, as illustrated in figure 2. 

The free-vibration problem is solved to obtain the natural frequencies 
^mn putting Pg(x,y,t) equal to zero in equation (Ale), combining equa- 

tions (Al) and (A5) , and solving the resulting simultaneous equations in u^, 
Vj^, and w^. The same frequency equation yielding the same set of natural 
frequencies O)^ is also obtained by combining equations (Al) and (A6) . 

Since both equations (A5) and (A6) satisfy the equations of motion and 
boundary condition of the fuselage structure, a complete solution for the 
response of the shell will involve the summation of both sets of solutions. 
Accordingly, the solutions for u, v, and w are taken as 



(X7a) 


(A7b) 


(A7c) 


where i = 1, 2. In equations (A7) , the subscript i = 1 corresponds to the 
set of solutions in equations (A5) and i = 2 corresponds to the set of solu- 
tions given in equations (A6) and 



(A8) 


(A9) 
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After combining equations (A7) and equations (Al) , the first two equations (Ala) 
and (Alb) identically vanish, and equation (Ale) , after simplifying and applying 
the orthogonal property of the normal modes, yields the equation of motion for 
the (mn i ) mode a s 




03 q . (t) 
mn^mni 


^mni 

M . 
mnx 


(AlO) 


since i = 1 and i = 2 yield the same natural frequencies The terms 

and are the generalized mass and force of the (mni) mode, respec- 

tively, given by 


and 

= 5J Pe(^'y't)W^^(x,Y)dx dy 


(All) 


(A12) 


The presence of light damping in the fuselage structure can be accounted for by 
adding a damping term in equation (AlO) so that the equation of motion of the 
(mni) mode of the damped structure becomes 


^^mn 


. (t) + 


03„ 


'^mni 




^mni('*^) 

%ni 


(A13) 


where Qj^ini denotes the dynamic magnification factor at resonance for the 
(mni) mo<ie. Equation (A13) is the same as equation (1) of the Analysis. 
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INTERIOR NOISE 

The internal acoustic field is defined by the wave equation 

and the boundary conditions. Assuming the process to be harmonic, the solution 
to the wave equation for acoustic pressure is 

p(x,(j),r,t) = jc^J^(ar) + C 2 Y^(ar^ 

X (cos Ticj) + C 3 sin r|<l5) (sin ijjx + cos i|^x) e^^^ (B2) 

where ^ constants of integration, and are hth 

order Bessel functions of the first and second kind, respectively, and 

a = kR[l - ^ 

r = (r/R) ) (B3) 

k = (ca/c) ^ 


and r, (p, and x denote the cylindrical coordinates along radial, circumfer- 
ential, and axial directions, and c is the speed of sound in air. 

For this problem open-end conditions are assumed for the fuselage, and 
hence the acoustic pressures are zero at x = 0 and x = L. Using these 
values in equation (B2) 


C4 = 0 

ilJ = ( 5 tt/l) 


(B4) 


For finite pressures at the center of the cylindrical section, C 2 must be 
zero. Hence, by summing up the contributions of all the modes, the interior 
pressure is obtained as 
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p(x,(f),r,t) = 2^ C^Ji^Cotr) j^in -^^(cos T]^ + C 3 sin (B5) 

?n 

where ^ - 1, 2, • . . , and d = 0/1/2, . . With the circiunferential 

coordinate 0 = (y/R) in equation (B5) , the total pressure in the interior of 
the fuselage can be expressed as 


p(x,y,r,t) 


?nj 


(B 6 ) 


where j = If 2 and 


= sin (Cttx/L) sin (ny/R) 
= sin (Cttx/L) cos (ny/R) 


(B7) 


When the process is harmonic, the pressure p is related to the particle dis- 
placement in the radial direction w^ by the relation 


w^ 


= ^1 /pQ03^j Op/3r) 


(B 8 ) 


where is the air density within the fuselage. This relation couples the 

'Structural and acoustic analyses. By equating the radial deflections of the 
shell and the particle displacement in the radial direction at r = R, it is 
possible to arrive at an expression for the interior noise pressure. Combining 
equation (B 6 ) and equation (B 8 ) the particle displacement in the radial direc- 
tion at r = R is 


^a ^) 


where ( ) ' = ) . 

w(x,y,t) in equation ( 8 ) , the values of the constants obtained 

by using the orthogonal property of the normal modes. It can be seen from 
equations (A7c) , (A 8 ) , and (A9) and equations (B 6 ) and (B7) that the space 

functions describing the shell radial deflections and the space functions 
describing the interior noise pressures at r = R are similar, that is. 




(B9) 


r . Po^ ^ 


Equating w^(x,y,R,t) given in equation (B9) to 




(BIO) 
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where mni = Cn j • Using this result and the orthogonality property of the 
space functions, the modal constants obtained as 


C 


mni 


PqOO^RPI 




(BID 


Substituting for from equation (Bll) in equation (B6) the noise pressure 

in the fuselage interior is obtained as 


|jn(a?)] (^12) 


In equation (B12) when ^ noise pressure goes to infinity. 'In an 

undamped acoustic field, subjected to finite excitations, pressure will go to 
infinity only at its natural frequencies. Hence, the roots of the equation 

J^(a) = 0 (B13) 

correspond to the natural frequencies of acoustic modes. If represents 

roots of equation (B13) , then the natural frequencies of acoustic modes are 
obtained using equations (B3) 


p(x,y,r,t) = PqOJ^rp ^ 


mni 




and 


= (mTTc/L) 


From equation (Bl2) the mean-square pressure is obtained as 


P 


2 


(x,y,r,03) 


2 2 4 
P PqR (jD 


mni 


2 4 

%ni^mn 



j^mni^^'^^ j^mnil J I (cross terms) (Bl5) 
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Space averaging the mean-square interior noise pressure, the cross terms in 
equation (B15) vanish due to the orthogonality of the modes, and hence the 
space -averaged mean-square pressure in the fuselage interior is obtained as 



for o)j^ < 0) where 


Y 


mni 


1 

27TRL 



mnx 


(x,y)dx dy 


(B17) 


Since a = kRj|l - (03j^/a))J ^ a 

case the equation (B16) becomes 


becomes imaginary when 


>0), in which 


2 , X 

P ((Ai) 


^2 2 2 4 
P PqR Cji) 


I 




'-mni 


2 4 
M .0) 
mni mni mn 


a 




Ymni 

¥ 


- Vl(l«l)Vl(l“D] 


(B18) 


> 0 ) where represents modified Bessel functions of the first kind. 

So far in the analysis, damping has not been considered for the acoustic 
modes, and hence the interior pressure will tend to infinity when J^id) ^ 0. 
Damping in the acoustic modes is introduced as explained in reference 10. The 
term Jj^(a) is expanded in a series around those points where Jj!^(a) = 0; 
that is. 


<<a) - + (o 


“ns' 




] 


(B19) 


Since ~ ®" 3 uation (B19) can be written as 


J (a) « (a - a )J (a ) 
n' ' ns n ns 


(B20) 
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where the series in equation (B19) is truncated after the second term. For 
small values of (a - ot^g) / 


- “L * 2ans(a - a^s) 


(B21) 


Substituting ^ 1 - = ‘is ~ 


xn 


equation (B21) it is seen that 


Jc^r2 


“ - “ns * 




2 a 


(B22) 


'ns 


Using the relationships in equations (B20) and (B22) , 


1 4 


j,4r4[-j..( )12 

1 - ( “ f 

L n '' ns'J 

\^mns/ 


(B23) 


where a ^ns* Inserting the acoustic dynamic magnification factor and 

allowing o) ^mns' equation (B2 3) can be written as 


1 


a 


2 



2 



(B24) 


Substituting the value of ( 1/a^ ) j^Jj^ (a)j given in equation (B24) into equa- 
tion (B16 ) , the space-averaged mean- square interior noise pressure in the 
vicinity of acoustic resonances is obtained as 


(W) 




H (0)/(jO ) 

' mn 


.,2 4 

. M .00 
mnx mnx mnx 


^mni I ^mni (^^ac) f 2 , ^ ^ i \ -r / 71 

-[Jn(a) - Jn_i(a)Jn+i(a)J 


1c4r4 


[jn(“ns>] 


(B25) 


for 00 ^ 00, 


mns 
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DETAILS OF THE FUSELAGES 


Fuselage 1 

Fuselage 2 
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TABLE II.- PROPERTIES OF MATERIALS 

Modulus of elasticity of skin, stringer, and frame, 

E,Ej.,Eg, Pa 0.71 X 10^^ 

Poisson's ratio, y 0.33 

Density of skin, stringer, and frame material, 

P,Ps,Pr, kg/m3 0.27 x lo^ 

Density of air, p^, kg/m^ 1.25 

Velocity of sound in air, c, m/sec 341.4 

Magnification factor for structure, Qj^ni 

Magnification factor for air, 25 
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Figure 3.- Typical circumferential and radial acoustic node 
patterns, n indicates number of circumferential full 
waves; s, number of radial waves. 


Frequency, Hz 


Theory 



Figure 4.- Comparison of predicted and measured natural frequencies 




(a) Fuselage 1. 


m 



(b) Fuselage 2 . 


Figure 5.- Modal frequencies of fuselage structure 
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Circumferential waves, n Circumferential waves, n 

(a) Fuselage 1. 

Figure 6.- Acoustic modal frequencies of fuselage interior space, m indicates number 
of longitudinal half waves; n, number of circumferential full waves; s, number of 
radial waves. 




Circumferential waves, n Circumferential waves, n 

(b) Fuselage 2. 


Figure 6.- Concluded. 
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Noise 

reduction, 

dB 



Frequency, f/f^ 

(a) Fuselage 1; exterior pressure distributed uniformly over 
the whole fuselage; f^ = 347 Hz. 


U) 

in 


Figure 8.- Variation of noise reduction with excitation frequency with and 

without sweep velocity. 




Frequency, f/f^ 

(b) Fuselage 2; exterior pressure distributed uniformly over 
the whole fuselage; f^^ = 1056 Hz, 


Figure 8.- Concluded, 







Frequency, Hz 

(b) Fuselage 2; exterior pressure on part of fuselage. 
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Figure 9.- Concluded. 





Sweep velocity, V , m/sec 

(b) Fuselage 2. 


Figure 10.- Concluded. 



Noise reduction, dB 



Frequency, Hz 

(b) Fuselage 2. 


Figure 11.- Concluded. 
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